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Lattice calculations of the QCD trace anomaly at temperatures T < 160 MeV have been shown 
to match hadron resonance gas model calculations, which include an exponentially rising hadron 
mass spectrum. In this paper we perform a more detailed comparison of the model calculations 
to lattice data that confirms the need for an exponentially increasing density of hadronic states. 
Also, we find that the lattice data is compatible with a hadron density of states that goes as 
p(m) ~ m~ a exp(m/Tff ) at large m with a > 5/2 (where Tb ~ 167 MeV). With this specific 
subleading contribution to the density of states, heavy resonances are most likely undergo 2-body 
decay (instead of multi-particle decay), which facilitates their inclusion into hadron transport codes. 
Moreover, estimates for the shear viscosity and the shear relaxation time coefficient of the hadron 
resonance model computed within the excluded volume approximation suggest that these transport 
coefficients are sensitive to the parameters that define the hadron mass spectrum. 
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I. INTRODUCTION 



Particle flow anisotropics at low transverse momentum produced in ultrarelativistic heavy ion collisions can be 
reasonably described [l|-y| using relativistic fluid dynamics with a very small shear viscosity to entropy ratio rj/s ~ 
l/(47r). This is as small as the uncertainty principle-based estimate derived by Danielewicz and Gyulassy nearly 
30 years ago [j] and also the more recent calculations [f| performed in strongly-coupled gauge theories dual to 
higher dimensional theories of gravity [f|. Beyond leading log perturbative QCD calculations that are applicable 
at temperatures T > m P i on give values for the ratio that are an order of magnitude larger than the bound [7J (for 
calculations based on parton transport see Q). Moreover, calculations performed using hadronic models at T ~ rn P i on 
also resulted in values for the ratio above the viscosity bound 

Due to the observation made in [l3[ that a small value of rj/s in QCD should occur in the transition region 
T ~ 150 — 200 MeV due to the rapid increase in the entropy density observed in lattice simulations [13, EH, the 
effects of an exponentially increasing density of hadronic states on several properties of hot hadronic matter were 
investigated using the hadron resonance gas model in Refs. [16-20]. It was shown in those studies that the addition 
of new hadronic states that follow an exponentially increasing hadron mass spectrum as proposed by Hagedorn (2lj 

e m/T H 

lim p(m) ry— (1) 

m— foe VTl * 

to the hadron resonance gas model led to a much better agreement to the lattice data computed around the pseudo- 
critical QCD critical temperature T c ~ 196 MeV inferred from Ref. (22[. Moreover, an estimate of rj/s at T ~ 190 
MeV computed using this model indicated that excited hadronic matter at those temperatures could become a nearly 
perfect fluid [23[ where rj/s approached l/(47r). In this case, the transition from viscous hydrodynamics to typical 
hadronic transport would be much smoother than expected. However, with the advent of the lattice calculations 
published in 
confirmed in 
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and the smaller critical region T ~ 155 MeV obtained in that study (which has been independently 
24j). a revision of the effects of heavy resonances on the hadron resonance gas model became necessary. 
Ref. [25j showed that a hadron resonance gas model containing only the known hadrons and resonances could only 
describe the lattice data up to T ~ 140 MeV while the inclusion of states with mass m > 2 GeV (which follow an 
exponential spectrum) could improve the match to the lattice data of [1 51 ] and provide a good description of lattice 
QCD thermodynamics up to T = 155 MeV. 

In this paper we present a more detailed comparison between the hadron resonance gas model calculations and the 
lattice data that not only provides strong evidence for the need of an exponentially increasing density of hadronic 
states with mass to > 2 GeV but also indicates that the density of states goes as p(m) ~ m~ a exp(m/T#)/ at large 
to with a > 5/2 (where Tjj ~ 167 MeV). As in (25j |. we estimate that the maximum temperature at which the 
hadron resonance gas model is applicable is ~ 155 MeV. A rough estimate of the shear viscosity computed within the 
excluded volume approximation for the hadron resonance model suggests that the shear viscosity to entropy ratio of 
hot hadronic matter is sensitive to the parameters that describe the hadron mass spectrum. We briefly comment also 
on the value of the shear relaxation time coefficient of hot hadronic matter. 



II. HADRON MASS SPECTRUM AND THE HADRON RESONANCE GAS MODEL 

In the statistical bootstrap model of hadrons pll . [26j . a hadron is considered to be a volume V (with typical 
length ~ 1 fm) composed of two or more freely roaming constituents and the hadron density of states is required 
to be consistent with the spectrum of constituents, which are themselves hadrons. This is the so-called "bootstrap 
condition" pionered in this context by Hagedorn in 1965 [211 ] . Frautschi (26| reformulated the bootstrap condition 
and wrote down an equation for the total density of states p(m) in the hadronic volume (for Boltzmann statistics) 
where different states of mass m s and energy Ei = \fpl + nif inside the box possess single particle density pi n (jni). 
The bootstrap condition 

lim p(m) => p in {m) (2) 

m— ^oo 

is exactly satisfied when limm-^oo p(m) ~ ce bm /m a with a > 5/2 26]. This exponentially rising mass spectrum is 
typical of a system consisting of string-like constituents such as a gas of free strings [27l - l29j or large N c glueballs [3(| ■ 
The spectrum of experimentally measured hadrons [3lT ] was found to be compatible with an exponential increase of 
the number of states up to ~ 1.7 GeV [HI, [33[ although the subleading power a cannot be reliably determined from 
such an analysis. 
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T H (GeV) 


A 


m (GeV) 


Q 


Pi 

P2 
P3 
P4 


0.252 
0.180 
0.175 
0.158 


2.84 (1/GeV) 
0.63 (GeV 3/2 ) 
0.37 (GeV 2 ) 
0.51 


0.5 
0.5 
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TABLE I. Parameters for the mass spectra shown in Eqs. ((8))- (|ll[) . 



Once the hadron mass spectrum is given, the thermodynamical quantities of the hadron resonance gas model (at 
zero chemical potential) in the total volume V and temperature T are fully determined by the partition function 
(assuming Boltzmann statistics) 
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which can be used to determine the usual thermodynamic functions. In fact, one finds 
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for the energy density, 
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for the entropy density, and 
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for the trace anomaly. The speed of sound can be found from the relation c 2 s 
different forms for the density of states 
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dpj de. In this paper we will discuss 4 
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where the parameters are shown in Table |U 

The parameters chosen for p\ and p 2 are the same used in [25j | while the parameters for the other p's were obtained 
from a fit to the lattice data of Ref. [15( . The density p% satisfies the asymptotic bootstrap condition exactly while 
pi (introduced in [HI]) and Hagedorn's p 2 satisfy the bootstrap condition within a power of m [2t| . The power law 
increase given by p4 (introduced by Shuryak in the early 70's [34]]) does not satisfies the bootstrap condition (it also 
does not lead to any singularities in the thermodynamics) but it provides a nice alternative to describe the rise of the 
hadronic mass spectrum. As discussed in (25|, the trace anomaly at temperatures around 160 MeV becomes sensitive 
to the heavy states in the spectrum with mass to > 2 GeV. A comparison between the different p's used here can be 
found in Fig. [1] where the integrand (in units of 1/GeV) in Eq. ([7]) is plotted as a function of to at T = 150 MeV. One 
can see that the integrand computed using the first three p's (black solid, blue dashed, red long dashed, respectively) 
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m (MeV) 

FIG. 1. Comparison between the integrand (in units of 1/GeV) in Eq. (0 computed using the different p's at T = 150 MeV as 
a function of m. The black line was computed using pi, the blue dashed line with p2, the red long dashed line with pz while 
the short dashed green line was obtained using p4. 



are very similar but they can be clearly distinguished from the result obtained using the power law in p± (small dashed 
green curve). 

Even though the upper limit of the mass integrals in Eqs. (j3][7]) is taken to infinity, the divergences implied by an 
exponentially rising spectrum do not appear in the calculations performed here because the limiting temperatures 
Th^s (i — 1,2,3) are above the largest temperature considered in this paper ~ 160 MeV. One may wonder if the 
approximations made in (j^JT]) (i.e., classical statistics and continuous mass spectrum) are at all justified. After all, 
we know that the measured hadroni c sp ectrum is of course discrete and that baryons and mesons obey different 
statistics. However, as pointed out in [25 1, the simplified formulas in Eqs. (|4][7]) provide an excellent description of the 
thermodynamic properties of the hadron resonance gas computed using the measured hadrons in the particle data 
book with the correct statistics in the temperature range T ~ 100 — 140 MeV if one imposes an upper cutoff for the 
mass integrals. In fact, for p\ and P2 the mass cutoff is 1.7 GeV and 1.9 GeV, respectively 25]. We have verified 
that p/T 4 computed using p$ with a mass cutoff of 1.9 GeV approaches the result obtained using p\ with a mass 
cutoff of 1.9 GeV. When T < 100 MeV, the discreteness of the hadron spectrum becomes relevant and the continuous 
approximation discussed here gives a poor description of the thermodynamic quantities of a hadron resonance gas. 
Therefore, in order to have a hadronic equation of state that is both valid at low temperatures (T < 100 MeV) and 
higher temperatures a hybrid model containing the measured hadron states plus a continuous Hagedorn spectrum 
above a certain mass cutoff is more appropriate [l6l - [2b1 | . 

A comparison between p(T)/T 4 , (e - 3p)/T 4 , and c 2 s {T) of the model defined by Eqs. dUtZJ) (for the four different 
hadron density of states) and the N t = 10 lattice data of Ref. [ll| can be found in Figs. [2KS] The black solid curves 
denote the result obtained by taking the mass integrals in Eqs. (J3KZJ) to infinity while the dashed blue curves were 
computed imposing an upper mass cutoff that varied for each p: for p\ the cutoff is 1.7 GeV while for pi and p3 the 
cutoff is 1.9 GeV. These cutoffs were determined by requiring that the trace anomaly computed in this continuous 
model matches the result (up to T ~ 140 MeV) obtained in a model where all the hadron resonances of the particle 
data book are included, as denned in (25[. 

The power law increase given by considerably simplifies the integrals in (j^JT]) and all of them can be done 
analytically. For instance, one finds in this case that c 2 s = l/(a + 4). However, note in Fig. [5] that the power law 
spectrum lacks the exponential growth necessary to describe the lattice data for temperatures above 140 MeV. This 
provides evidence that the thermodynamic quantities of QCD computed on the lattice can be understood in terms of 
a simple hadron resonance gas model with an exponentially rising density of states at temperatures T ~ 100 — 155 
MeV. This comparison to lattice data cannot pin down the exact subleading power of m in the hadron density of 
states. However, it is important to note that the specific value of this power has some interesting consequences, as we 
shall elaborate below. 

The subleading contribution ~ mT a to the density of states at large m, according to Frautschi's seminal paper 
[2fj| . determines the decay properties of a heavy resonance. For instance, when a > 5/2 (which is the case of p%), a 
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heavy resonance decays (in the first generation of its decay chain) into a heavy secondary particle that carries almost 
all the available mass plus one (with 69% probability) or two light hadrons (with 24% probability). This should be 
contrasted with Hagedorn's original mass spectrum in Eq. ([9]) for which the statistically favored process involves a 
heavy resonance of mass m decaying into a number n ~ In m of secondary particles, each of similar mass (26|. 

It is well known that hadrons interact with each other in a variety of different channels, some of them give rise 
to repulsive interactions while others represent attractive interactions. In Ref. (35j the pressure of an interacting gas 
of pions was calculated within the virial expansion (using experimentally determined phase shifts) and it was shown 
that the thermodynamic quantities of this interacting system nearly coincides with those of a free gas of pions and 
p mesons. In this case, there is an approximate cancellation between the attractive and repulsive S-wave channels 
which effectively enhances the P-wave contribution from the p resonance [36[. As more hadronic species are included, 
it is not at all guaranteed that the standard assumption behind hadron resonance models, i.e., that the interacting 
hadronic system can be described by a free gas of the original hadrons and their resonances, is applicable. In general, 
the inclusion of resonances represents the contribution from the attractive channels while repulsive interactions can 
be modeled using simple excluded volume corrections to the thermodynamics [37i - 439j . The suggestion, obtained 
from a comparison to lattice data, that the complicated interactions among hadrons that enter in the calculations of 
QCD thermodynamics at temperatures of the order of the pion mass can be effectively modeled by a simple, non- 
interacting gas of hadrons and resonances in accordance with the bootstrap model [2l[ is therefore quite unexpected 
and remarkable. 

Given the known uncertainties in lattice calculations at low temperatures, the conclusion made above regarding the 
applicability of the hadron resonance gas should be taken with great care. If results obtained with finer lattices confirm 
this picture, this would provide very strong evidence for the validity of Hagedorn's bootstrap hypothesis. We note in 
passing that recent lattice calculations of the thermodynamical properties of SU (3) pure glue have shown evidence for 
the presence of an exponentially rising glueball mass spectrum below the deconfinement critical temperature [4(J 5l[ . 



III. EXCLUDED VOLUME CORRECTIONS TO THE HADRON RESONANCE GAS MODEL 

As mentioned in the previous section, in general one should expect that there are repulsive interactions among 
hadrons and that a simple way to take that into account in the hadron resonance model is via the excluded volume 
corrections (37- 39] . In this case, the partition function in (j3|) becomes 

(\ N 
V-jlv] (12) 
i=i / 

where Vj denotes the excluded volume by the jth hadron. We shall assume for simplicity that the volume excluded 
by each hadron is a constant that is basically the same for all hadrons, i.e., Vj — v. This parameter can be written 
in terms of an effective hard-core volume, v = 4 • 4-7rr 3 /3, where r is the effective core radius. The excluded volume 
pressure is determined by the equation 

P^l=n(T)eM-vp v (T)/T) (13) 

where n(T) — p(T)/T is the total particle density computed without volume corrections (39l |. The equation above 
can be solved analytically in terms of the Lambert W function [42[ and it reads 

E^l = W(vn(T)). (14) 

The other thermodynamic quantities, e v (T), s v (T), and n v (T) can be obtained through the pressure using the standard 
thermodynamic identities. In the limit where v — > one recovers the formulas in (UK]). 

Given that the free hadronic gas of the previous section provided a good description of the data, one should expect 
that the volume corrections should be minimal in this case. In fact, one can again use the lattice data and the different 
hadron mass spectra discussed before to show that the excluded radius cannot be larger than 0.2 fm for pi (i = 1, 2, 3). 
Pi is not considered in this case since it fails to describe the data for T > 140 MeV. We show a comparison between 
the lattice data and the p% model curves for r = 0.2 fm in Fig. [6] The small excluded volume shifts the curves slightly 
downwards, which makes them get closer to the lattice data at T = 160 MeV. Similar results hold for the other 
exponentially increasing spectra considered before. For larger excluded volumes the hadron resonance gas curves 
start to deviate from the lattice at lower temperatures and this is why we here take r < 0.2 fm. This analysis shows 



6 



e - 3 p 




100 110 120 130 140 150 160 170 
T (MeV) 



1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



: — M->oo 


A 


1 --M^1.7 GeV 


>• ✓ 




Pi : 



100 110 120 130 140 150 160 170 
T (MeV) 




0.00 



100 110 



120 130 140 
T (MeV) 



150 



FIG. 2. Trace anomaly, pressure, and speed of sound squared for the hadron resonance model with density of states pi. The 
black solid curves denote the result obtained by taking the mass integrals in Eqs. (|4I7|) to infinity while the dashed blue curves 
were computed imposing an upper mass cutoff of 1.7 GeV. The data points correspond to the Nt = 10 lattice data published 
in Ref. [IS] (obtained from table 5 in that paper). 
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FIG. 3. Trace anomaly, pressure, and speed of sound squared for the hadron resonance model with density of states p2- The 
black solid curves denote the result obtained by taking the mass integrals in Eqs. (|4I7|) to infinity while the dashed blue curves 
were computed imposing an upper mass cutoff of 1.9 GeV. The data points correspond to the Nt = 10 lattice data published 
in Ref. [IS] (obtained from table 5 in that paper). 
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FIG. 4. Trace anomaly, pressure, and speed of sound squared for the hadron resonance model with density of states ps. The 
black solid curves denote the result obtained by taking the mass integrals in Eqs. (|4I7|) to infinity while the dashed blue curves 
were computed imposing an upper mass cutoff of 1.9 GeV. The data points correspond to the Nt = 10 lattice data published 
in Ref. [IS] (obtained from table 5 in that paper). 
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FIG. 5. Trace anomaly, pressure, and speed of sound squared for the hadron resonance model with density of states pi. The 
black solid curves denote the result obtained by taking the mass integrals in Eqs. (07} to infinity. The data points correspond 
to the Nt = 10 lattice data published in Ref. [15| (obtained from table 5 in that paper) . 
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that volume corrections do not play a significant role in the description of lattice data (at least for the mass spectra 
parameters determined in the previous section). Of course, given the known uncertainties in lattice calculations at low 
temperatures, another possibility would be to define the parameters in a way that the model including the continuous 
spectrum fits the lattice data only at higher temperatures around T ~ 150 MeV, as it was done in [17|. This will be 
discussed further in the next section. 



, 1/2 

V = 64 r 2 



IV. CALCULATION OF THE SHEAR VISCOSITY TO ENTROPY DENSITY RATIO AND THE 
SHEAR RELAXATION TIME COEFFICIENT OF THE HADRON RESONANCE GAS MODEL WITH 

EXCLUDED VOLUME CORRECTIONS 

The computation of the transport properties of a hadronic mixture is not an easy task. There have been several 
studies on this subject in the last few years 0412) G3, EH EH- In order to find at least an estimate of the order of 
magnitude of the rj/s ratio of hot hadronic matter at T ~ 160 MeV, we follow the approximations made by Ref. @ 
where the shear viscosity of a multi-component gas of hadrons and resonances in the excluded volume approximation 
(described in the previous section) is given by 

/ rp \ 1/2 rp „oo /m\ 

{-) dm P (m)m^K 5/2 (-) . (15) 

Currently, it is not known how to compute the contribution to the shear viscosity from heavily massive and highly 
unstable resonances that cannot rigorously be described using the Boltzmann equation. These states contribute 
significantly to the thermodynamic properties of the matter at high tempetatures (as shown in the previous sections) 
and, due to their rapid decay, it is natural to assume that their presence will affect the mean free paths of the other 
hadrons. In [ijj , it was assumed that the mean free path of these resonances with m > 2 GeV equals their inverse 
decay width. Obviously, further studies have to be carried out to properly include the effects of Hagedorn states on 
transport coefficients of hot hadronic matter. While the formula in Eq. (j!5[) may only provide an estimate of the shear 
viscosity of an interacting hadron gas, the temperature behavior of r\ computed with this approximation Q follows 
the estimates made using other methods and thus we shall proceed using this formula. 

In Eq. (|15p. the dependence on the excluded radius only appears via the 1/r 2 factor. However, when computing 
rj/s, the entropy density should be the one determined within the same approximation, i.e, s v (T). Therefore, even 
in this approximation rj/s possesses a nontrivial dependence on the hadron cross section ~ r 2 . This is shown in 
Fig. [7] where the rj/s ratio is computed for different r's using p% (|T0|) and T = 155 MeV. Note that for the highest 
temperature considered T = 155 MeV, rj/s depends very weakly on r when the excluded radius is larger than 0.2 fm 
(using this specific set of parameters that define the mass spectrum). 

In Fig. [5] we show the temperature dependence of rj/s computed using the density of states p% in Eq. (|TU1) and 
r = 0.2 fm. As discussed in the previous section, for the specific choice of parameters that define ps, calculations 
of thermodynamical quantities performed with an excluded radius of r — 0.2 fm can describe the lattice data in the 
entire temperature range, T = 100 — 160 MeV. This value of the excluded radius is smaller than other estimates 
@, which however were not constrained by fitting the lattice data. Perhaps it would be more physical to consider a 
model where the excluded radius increases with the mass of the hadron but for simplicity's sake in the current we 
will limit ourselves to a constant excluded radius hoping that we are correct within an order of magnitude. In Fig. [FJ 
we see that the rj/s ratio remains an order of magnitude above the viscosity lower bound up to T = 155 MeV (this 
remains the case when other expressions for the density of states mentioned in the previous sections are used). The 
entropy density computed in this case matches the lattice data well, as it can be inferred from the other quantities 
shown in Fig. ©. Therefore, this simple hadron resonance gas model with constant excluded volume corrections is 
able to describe the thermodynamic quantities computed by lattice data below T = 160 MeV and the corresponding 
r\ is computed self-consistently within the same framework. 

In viscous hydrodynamic calculations of the QGP time evolution (TMl], there is at least another transport coefficient 
that must be included in the fluid equations, the shear viscosity relaxation time T n , which enters in second order 
viscous hydrodynamic calculations. In fact, in relativistic fluids causality is intimately connected to stability [45|, Efi| 
and Israel and Stewart [47| were among the first to understand that the characteristic times within which fluid 
dynamical dissipative currents relax towards their asymptotic Navier-Stokes values cannot be arbitrarily small. Using 
the Boltzmann equation, it is possible to show [47H50| that in relativistic gases t v is of the order of the microscopic 
collision time. A detailed linear stability analysis made in Ref. 51] showed that stability and causality require that r v 
in any viscous relativistic fluid cannot be arbitrarily small. In fact, this transport coefficient must obey the following 
inequality [5l| 

4 V 1 

r ^3lr(T^)- (16) 
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FIG. 6. Trace anomaly, pressure, and speed of sound squared for the hadron resonance model with density of states p3 and an 
excluded volume radius of r = 0.2 fm. The dashed blue curves denote the excluded volume results while the solid black curves 
show the corresponding quantities without excluded volume corrections. The data points correspond to the Nt = 10 lattice 
data published in Ref. [Tg] (obtained from table 5 in that paper). 
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FIG. 7. The ratio rj/s (dashed blue line) as a function of the excluded volume radius computed via (|15p using the density of 
states pz in Eq. (|10|) and T = 155 MeV. The black line denotes the viscosity lower bound r\js — l/(47r). 
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FIG. 8. The ratio r//s (dashed blue line) computed within the excluded volume approximation defined via Eqs. (|13fl and (1151) 
using the density of states ps in Eq. (|10[) and r = 0.2 fm. The entropy density computed in this model fits the lattice data [15|] 
in the temperature range shown. The black line denotes the viscosity lower bound rj/s = l/(47r). 
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FIG. 9. The lowest relaxation time coefficient tv computed using the inequality in Eq. (|16[l for the hadron resonance gas model 
defined with density of states p3 and rj/s and c? s computed within the excluded volume approximation with excluded volume 
radius r — 0.2 fm (dashed blue curve). The black line denotes the lowest value for r n in a conformal "nearly perfect" fluid 
where (? s = 1/3 and rj/s = l/(4m). 



While we shall not compute tv for the hadron resonance gas model considered here, we find it instructive to consider 
the smallest tv implied by the inequality above since it provides an estimate for the value of this parameter that can 
be used in hydrodynamic simulations. In Fig. @ we show (dashed blue line) the lowest value for tv computed using 
the rj/s in Fig. [8] that fulfills the stability and causality criteria. We also show in the same plot the lowest value for 
7v in a generic conformal "nearly perfect" fluid where cj, = 1/3 and rj/s = l/(47r) for reference. Note that the lowest 
value of tv computed in the hadron resonance gas model employed here remains well above the lowest value for a 
nearly perfect fluid up to T — 160 MeV. 

Throughout our calculations thus far we have the underlying assumption that the hadron gas model must fit the 
entire lower temperature region T ~ 100 MeV of the lattice data. However, it is interesting to consider the possibility 
of just fitting the lattice data at higher temperatures T ~ 150 MeV to see how that affects the rj/s calculation. This 
can be accomplished by changing some of the parameters of the hadron mass spectrum. Before doing so, we increase 
r = 0.5 fm Q. Then we refitted the parameters ^3 and rrios in (I10p in order to fit only the high temperature region 
of the lattice trace anomaly. This is shown in Fig. (1101) . The corresponding rj/s and shear relaxation time computed 
with this setup are shown in Figs. (fTTj) and ([12")) . Note that with this new set of parameters rj/s drops down to l/(47r) 
around T = 160 MeV and r w also decreases significantly in that region. The decrease of the transport coefficients 
is mainly driven by the larger entropy that results from this new fit that is only constrained by the higher lattice 
temperatures. This emphasizes the sensitivity of such calculations to the thermodynamic properties of the matter. 



V. CONCLUSIONS 



In this paper we performed a detailed comparison of hadron resonance gas model calculations to recent lattice 
data (TBI that confirmed the need for an exponentially increasing density of hadronic states. Also, we showed that 
the hadron mass spectrum, extracted from a comparison to lattice data, is compatible with p(m) ~ mT a exp(r7i/Tff ) 
at large m with a > 5/2 (where Tjj ~ 167 MeV). With this specific m~ a (with a > 5/2) factor in the density of 
states, heavy resonances most likely undergo 2-body decay (instead of multi-particle decay) (26[, which facilitates 
their inclusion into hadron transport codes. Moreover, we have computed the shear viscosity to entropy density ratio 
of this system within the excluded volume approximation and the results suggest that rj/s of hot hadronic matter 
is very sensitive to the temperature dependence of the thermodynamic quantities. Using this calculation of rj/s, we 
were able to compute the lowest value for the shear relaxation time coefficient used in second order hydrodynamic 
calculations that respects the criteria of causality and stability of a relativistic viscous fluid. 

Previous estimates for the rj/s ratio in a hadronic gas (l7| had concluded that hadronic matter at temperatures 



14 



e - 3 p 




100 110 120 130 140 150 160 170 
T (MeV) 

FIG. 10. Comparison between the model's trace anomaly computed different excluded volumes and the Nt = 10 lattice data 
published in Ref. [l5[ (obtained from table 5 in that paper). The solid black curve was computed using p3 with the parameters 
defined in the table [I] The short-dashed blue curve is computed using the same p3 but with an excluded volume of r = 0.2 
fm. The long-dashed red curve was computed using p3 but with A3 = 3.7 GeV 2 and mo3 = 1 GeV together with an excluded 
volume of r = 0.5 fm. 
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FIG. 11. The ratio rj/s computed within the excluded volume approximation. The short-dashed blue curve is computed using 
P3 with the parameters defined in the table [T] but with an excluded volume of r = 0.2 fm. The long-dashed red curve was 
computed using P3 but with A3 = 3.7 GeV 2 and mo3 = 1 GeV together with an excluded volume of r = 0.5 fm. 



T ~ 190 MeV behaved as a nearly perfect fluid. This is not at odds with the findings presented in this paper 
as we explained in the previous section. In fact, the curve shown in Fig. (jHJl approaches l/(47r) when continued 
to temperatures ~ 190 MeV. Moreover, if only the high temperature region of the lattice data is fitted, then the 
transport coefficients calculated here decrease significantly (e.g., 77/s ~ 1/(4tt) near T — 160 MeV). This highlights the 
importance of knowning the correct temperature dependence of the thermodynamic quantities of QCD at temperatures 
~ 100 — 160 MeV since it may play a important role in the calculation of transport coefficients. 

The key difference between these studies is the lattice data used as a reference for the hadron gas calculations. Ref. 
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FIG. 12. The lowest relaxation time coefficient T n computed using the inequality in Eq. (|16[) for the hadron resonance gas 
model. The short-dashed blue curve is computed using ps with the parameters defined in the table U but with an excluded 
volume of r = 0.2 fm. The long-dashed red curve was computed using p3 but with A3 — 3.7 GeV 2 and mo3 = 1 GeV together 
with an excluded volume of r = 0.5 fm. 



17] used the most recent lattice data, at the time [22] which indicated a phase transition pseudo-critical temperature 
T c ~ 196 MeV. The much lower value for this pseudo-critical temperature found in fl5j severely reduced the value of 
the maximum temperature at which the hadron resonance gas is still applicable since "T c " decreased from 190 MeV to 
160 MeV. Given that this low pseudo-critical temperature has already been independently confirmed by other lattice 
groups 24| , if there is no change in the low temperature behavior of the thermodynamic quantities as determined by 
lattice, the analysis performed in this paper indicates that the hot hadronic matter formed in ultrarelativistic heavy 
ion collisions is far from being a nearly perfect fluid. However, this should be taken with a grain of salt given the 
sensitivity mentioned above within the transport coefficients. 

Since (e — 3p) /T 4 in the lattice data [l5| continues to increase until it reaches a turning point around T ~ 200 MeV, 
one may wonder if it is possible to devise a model that reduces to the Hagedorn resonance gas discussed here at low 
temperatures while also incorporates the correct degrees of freedom in the crossover region between T ~ 160 — 200 
MeV. Perhaps such an effective model can be constructed by taking into account the Polyakov loop [52[ . 
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Pesquisa do Estado de Sao Paulo (FAPESP). J. Noronha thanks Conselho Nacional de Desenvolvimento Cientifico e 
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